Electrostatic model of atomic ordering in complex perovskite alloys 
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We present a simple ionic model which successfully reproduces the various types of compositional 
long-range order observed in a large class of complex insulating perovskite alloys. The model as- 
sumes that the driving mechanism responsible for the ordering is simply the electrostatic interaction 
between the different ionic species. A possible new explanation for the anomalous long-range order 
observed in some Pb relaxor alloys, involving the proposed existence of a small amount of Pb"''* on 
the B sublattice, is suggested by an analysis of the model. 

PACS numbers: 64.60. Cn, 81.30.-t, 77.84.Dy, 77.84.-s 
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Complex insulating perovskite alloys are of great in- 
terest for their actual or potential uses mainly because 
of their exceptional dielectric and piezoelectric proper- 
ties [ij. Examples are the so-called "super Q" mixed 
metal perovskites, such as Ba(Zni/3Ta2/3)03 (BZT) 
or Ba(Mgi/3Nb2/3)03-BaZr03 (BMN-BZ), which can 
be used in high frequency applications because of 
their ultra- low losses at microwave frequencies [^D- 
Other examples are the "new" relaxor ferroelectric al- 
loys, such as Pb(Mgi/3Nb2/3)03-PbTi03 (PMN-PT) 
or Pb(Zni/3Nb2/3)03-PbTi03 (PZN-PT), which exhibit 
extraordinarily high values of the piezoelectric constants 
Q and thus promise to improve the resolution and range 
of ultrasound and sonar listening devices 

An aspect that is not well understood about these com- 
plex materials is the compositional atomic ordering that 
occurs in some compounds and not in others, and the 
different kinds of long-range order (LRO) that can oc- 
cur. For example, increasing progressively the bulk com- 
position a; in (1 — x) Ba(Mgi/3Nb2/3)03 + x BaZr03 
(BMN-BZ) leads to the following types of order for a 
final sintering temperature of 1640°C the so-called 
"1:2" ordering along the [111] direction for very small 
values of x (i.e., lower than 5%); then "1:1" ordering 
along [111] for compositions ranging from 5 to 25%; and 
finally, disordering for larger values of x. 

The nature and strength of the compositional order 
may have crucial consequences for the desired proper- 
ties of the material. For example, previous research has 
shown that the microwave loss properties of perovskite 
ceramics are very sensitive to the B-site cation order with 
the Q value improving with increasing degree of order 
Similarly, compositional fluctuations are believed 
to play a central role in the "relaxor" behavior in Pb- 
based perovskite alloys ]0] . 

Previous pioneering theoretical work ]^,^ has been 
partially successful in describing the ordering in some 
cases. In such models, the Coulomb interactions are in- 
cluded indirectly via a cluster expansion of finite range. 
While there is some justification for such an approach 
[pO|, the ordering sequences that occur in complex ma- 



terials such as BMN-BZ did not emerge naturally, and 
thus the driving mechanism responsible for the ordering 
has remained somewhat obscure. 

Motivated to resolve these problems, we decided to in- 
vestigate a simple model that explicitly includes the long- 
range Coulomb interactions between ions. In this model, 
the energy is taken to be proportional to the electrostatic 
energy of an ideal system of ionic charges. We initially 
expected that this very simple model might need to be 
augmented with corrective covalent terms before it could 
provide a useful account of the different orderings seen 
in real materials. To our surprise, the minimal model by 
itself is remarkably successful. The purpose of this let- 
ter, then, is to show that this very simple model provides 
a systematic understanding of the complicated ordering 
behavior of complex perovskite alloys, thereby providing 
strong evidence that the Coulomb interaction between 
ions is the dominant factor in determining such ordering. 

Our model assumes that the minimization of the elec- 
trostatic interactions between the different ions in the 
compounds leads to the LRO seen experimentally. The 
interaction energy in ABO3 perovskite compounds is 
then expressed as 
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where R;,- is the position in the ideal cubic structure of 
the atom on site r of cell / (t={A,B,Oi, 02,03}), and e 
is a dielectric constant providing some screening effects. 
While we find that this model can be successfully ap- 
plied to other cases (see below), we first limit ourselves 
to divalent-A compounds having a common A atom on 
all A sites, and compositional freedom on the B site. In 
this case, we have Qi^a = Qa = +2 and Qi^o = Qo = ^2 
independent of I. Focusing on the B sites on which the 
alloying occurs, we can decompose Qi,b as 



(2) 



where the average valence qb is equal to -1-4 while Aqi de- 
pends on I. In order to preserve overall charge neutrality, 
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the average Ag of Aqi over all cells must vanish. 

Substituting expression (2) for Qi^b into Eq. (1), we 
can write E = Eq + Ei + E2, where the subscript refers 
to the number of powers of Aq appearing in that term. 
Then Eq is just a constant; it is the energy, within 
this model, of an ideal divalent-A and tetravalent-B per- 
ovskite compound. Ei can be written as 
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but noting that the term in brackets does not actually 
depend on / and using Ag = 0, it follows that Ei vanishes. 
Thus the only term depending on the atomic order is 
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In fact, using R; — la where a is the lattice constant. 
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Eq. (5) demonstrates that our model can be reduced to 
focus only on the alloying sublattice (e.g., the B sub- 
lattice in the present case). Moreover, the model is 
parameter-free in the sense that e^/2ea defines a char- 
acteristic energy that will set the temperature scale on 
which ordering will occur. 

In the following, we apply our model to various cases of 
divalent-A compounds with B-site alloying. We thus in- 
troduce a B-site-only notation. For example, IVa;IV'i_a: 
will denote a homovalent binary alloy having tetrava- 
lent B atoms, e.g., Pb(Zr2;Tii_j;)03. In the same man- 
ner, n(i_a;)/3lVa;V2(i-2;)/3 will indicate a heterovalent 



ternary, such as (l — x) Ba(Mgi/3Nb2/3)03 + x BaZr03, 
which is the result of mixing a II1/3V2/3 heterovalent 
binary (Ba(Mgi/3Nb2/3)03) with a tetravalent-B bulk 
(BaZr03). 

In view of Eqs. (2) and (5), two points are imme- 
diately obvious, (i) The present model does not allow 
IVa;IV'i_2; homovalent binary alloys to order, since Aqi 
is identically zero in all cells I. This is in agreement 
with the experimental absence of LRO in homovalent al- 
loys ll^. (ii) II1/2VI1/2 and III1/2V1/2 heterovalent bi- 
naries should behave much the same, since Aqi consists 
in both cases of equal populations of charge ±n, where 
n = 2 and 1 for II1/2VI1/2 and III1/2V1/2 alloys, respec- 
tively. Thus, these two cases really only differ in the 
magnitude of the effective energy scale factor n^e^/2ea. 
By computing the energy of a large variety of candi- 
date structures and by performing Monte-Carlo simula- 
tions for small temperatures (see below), we determined 
that the ground-state (zero-temperature) structure pre- 
dicted by the present model in these two cases is rocksalt, 
i.e., an alternation of ±rt charged planes along the [111] 
direction. This is precisely the LRO observed experi- 
mentally in Pb(Sci/2Tai/2)03 (7|jl|, Pb(Ini/2Nbi/2)03 



y, Pb(Ybi/2Tai/2)03_Jl|, Pb(Mgi/2Wi/2)03 y, 
and Pb(Coi/2Wi/2)03 [y7^ 

Similar calculations show that the present model 
also reproduces the ground-state 1:2 structure ob- 
served in many II1/3V2/3 heterovalent binaries, such 



as Ba(Mgi/3Nb2/3)0, 



(BZT) 

PI, Ba(Sri/3Ta2/3, 
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(BMN) 
Ba(Mg 
(BST) 



Ba(Zni/3Ta2/3 



Ta2/3)03 



)03 

(BMT) 



Ba(Zni/3Nb2/3)03 

(BZN) |lj], Ba(Nii/3Nb2/3)03 (BNN) Q, and 
Ca(Cai/3Ta2/3)03 (CCT) ||. That is, we do indeed 
find that the structure minimizing the electrostatic en- 
ergy consists of one (111) layer of —2 relative charges 
(e.g., Zn-derived in BZN) alternating with two (111) lay- 
ers of -1-1 relative charges (e.g., Nb-derived in BZN) along 
the [111] axis. 

To extract further consequences from the model, it is 
necessary to solve for its thermodynamic behavior as a 
function of temperature. This is done using Metropolis 
Monte Carlo simulations (20| in which the (electrostatic) 
energy E2 is given by Eq. (5). To be definite, we as- 
sign a=7.7 a.u. and e=10. We normally limit ourselves 
to 6 X 6 X 6 supercells since we checked that increasing the 
size of the calculations up to 12x12x12 supercells does 
not change the results appreciably. We typically use a 
very large number of trial moves (up to 3 x 10^ in tran- 
sition regions) in order to obtain good statistics for each 
temperature. Since this corresponds to a large number 
(^35,000) of sweeps, the initial configuration is unim- 
portant; it can either be chosen at random, or set to a 
given atomic order such as 1:2 or 1:1. To find the ground 
state, we first perform calculations at high temperature 
(e.g., 4,000K), and then very slowly decrease the temper- 
ature until reaching a structure for which no moves are 
accepted by the Monte-Carlo program. 

Once the thermodynamic state corresponding to a 
given temperature has been reached, the search for LRO 
is done by calculating, from the Monte-Carlo outputs, the 
Fourier transform of the charge-charge correlation func- 
tion, 



7(k) ^ a^Aqi Aqi+v cxp(-ik • 1') 



(6) 



Here a is a normalization factor, the sum runs over the 
charges in the B sublattice, and k is the wavevector in 
the Brillouin zone of the unit cubic cell. It can easily be 
demonstrated that 77(k) is directly proportional to the 
ensemble average of the square of the Fourier transform 
of the charge distribution. Thus, a large value of rj at 
k — 27r(i,i,i) corresponds to strong 1:1 LRO along 
the [111] direction, while a peak in 77 at k = 27r(i, i, i) 
indicates 1:2 order along the [111] direction. 

Figures l(a-c) show the predicted behavior of 
77(k) vs. tetravalent composition x at T^IOOOK 
for III(i_:r)/2lVa;V(i_2.)/2, l\i-x) /2^^ x^\i-x) /2, and 
II(i_a;)/3lVa;V2(i_a;)/3 hctcrovalent ternaries, respec- 
tively. We find that progressively increasing x in 
in(i_a;)/2lV:EV(i_a;)/2 and ll{i^x)/2^^x^l(i~x)/2 leads to 
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FIG. 1. Monte Carlo simulations of the long-range order 
parameter ry(k) vs. tetravalent atomic composition x for 
the model of Eq. (5), with 6x6x6 supercells, a=7.7 a.u., 
e=10, and T^IOOOK. Solid and dashed lines refer to 1:1 order 
[k = 27r(i, i, i)] and 1:2 order [k = 27r(i, §, i)], respectively, 
(a) III(i_^)/2lVa;V(i_^)/2 ternaries, (b) II(i_^)/2lV^VI(i_^)/2 
ternaries, (c) ll(i-x)/3'^^ x'V2(i-x)/3 ternaries. 

a continuous transition from rocksalt-type order to a dis- 
ordered state. This is consistent with the experimental 
findings in (l-x) Pb(Sci/2Tai/2)03 + x PbTiOg |l9|,0. 
We also predict that the 1:1 rocksalt order can survive in 
II(i_2:)/2lVa;VI(i_j.)/2 alloys up to rather large x. This 
compositional difference between ni(-]^_^)/2lVa;V(i_2;)/2 
and /2IV2: VI(i_a;) /2 ternaries is due to the fact that 

the ground state energy of the latter is four times deeper 
than that of the former (for a given e). Thus, quite a 
large number of tetravalent atoms is required to induce 
a disordered state in the ll{i-x) /2i^ x^\i-x) /2 ternaries. 
We varied the temperature for the ll\i-x) /2i^ (i-x) /2 
and n(i_2,)/2lVxVI(i_2.)/2 ternaries but found no quali- 
tative changes (e.g., new phases), although naturally the 
transition from the ordered to the disordered phase is 
shifted to lower x with increasing T ||2^ . 

Figure 1(c) demonstrates that the atomic ordering pre- 
dicted in II(i_^)/3lV2;V2(i-a;)/3 ternaries is even richer 
than in the previous two cases. We find 1:2 order along 
[111] for small x, then 1:1 rocksalt order for interme- 
diate X, and finally a disordered state for even larger 
X. These predictions are in remarkable agreement with 
recent experimental observations of precisely this or- 
dering sequence as a function of x in the microwave 
ceramics (1 - x) Ba(Mgi/3Nb2/3)03 + x BaZr03 [|, 



(1-a;) Ba(Mgi/3Ta2/3)03 + x BaZr03 ||, and (l-x) 
Ba(Zni/3Ta2/3)03 + x BaZr03 We find numeri- 

cally that the critical value of x at which the crossover 
occurs from 1:2 to 1:1 order depends on the tempera- 
ture used in the Monte-Carlo simulation, as well as on 
the value of the dielectric constant in Eq. (5). Increas- 
ing the temperature or increasing the dielectric constant 
leads to a decrease of this crossing composition. For ex- 
ample, keeping the lattice constant equal to 7.7 a.u. and 
using a dielectric constant of 15 yields a critical com- 
position Xc ~ 5% for T=1000K. (A further increase in 
temperature or dielectric constant leads to alloy disorder 
even at small x composition.) Moreover, our calculations 
also confirm a suggestion based on experimental observa- 
tions jlJI^ for the structure responsible for the 1:1 order. 
Specifically, when we observe 1:1 order for x < 25%, we 
find a rocksalt structure in which one sublattice is almost 
entirely occupied by pentavalent atoms, while the sec- 
ond sublattice is basically composed of the divalent and 
tetravalent atoms plus the remaining pentavalent atoms. 
The figure also indicates that the strongest 1:1 order oc- 
curs for a bulk composition of 25%, which is consistent 
with the experimental findings that the l:l-order X-ray 
reflections seems to reach a maximum for x = 0.25 in 
(1 - x) Ba(Zni/3Ta2/3)03 + x BaZr03 (BZT-BZ) [||. 

On the other hand, there is a noteworthy exception 
to the ordering sequence shown in Fig. 1(c). Namely, 
in II(i-a;)/3lV2:V2(i_2;)/3 ternaries for which the diva- 
lent A atom is Pb, the 1:2 ordering is not observed 
as X — > 0. Instead, these compounds exhibit weak 
1:1 X-ray reflections all the way down to a; = 0, 
as observed in Pb(Mgi/3Nb2/3)03 (PMN) and 
Pb(Mgi/3Ta2/3)03 (PMT) As this appears to be 

at variance with our model, some explanation is clearly 
required. First, we note that our model predicts that the 
1:2 and 1:1 ordering types are very close in free energy 
in n(]^_^)/3lVa;V2(i_2;)/3 ternaries at small x. It is thus 
possible that in TL1/3V2/3 binaries for which the A sub- 
lattice is occupied by Pb, some physical mechanism may 
occur that is responsible for triggering the 1:1 order in 
place of the 1:2 type. For example, a possible mechanism 
might be the covalency induced by the existence of short 
Pb-0 bonds in the lead compounds |^^. The difference 
in covalency between Pb compounds and Ba compounds 
has been previously proposed to explain the difference in 
the ferroelectric behavior of PbTi03 and BaTi03 [|0). 

However, we wish to propose another possible mech- 
anism for the existence of weak 1:1 order in PMN and 
PMT, related to the multivalent nature of the Pb atom. 
As is well known, Pb can be either divalent or tetrava- 
lent, in view of its s^p^ valence electronic configuration. 
Indeed, as much as 5% of Pb+'* has been experimentally 
reported in PMN These tetravalent Pb ions were 

initially thought to be sitting on the A sublattice of the 
perovskite alloy |3^|. However, very recent experiments 
on Pb(Sci/2Tai/2)03 (PST) films indicate the presence 
of some Pb^^ on the B sublattice of the perovskite struc- 
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ture Thus, we would like to raise the possibility 

that the Pb"*"^ atoms in PMN sit on the B sublattice. In 
this way, the application of Eq. (2) to PMN must rec- 
ognize the existence of three species on the B sublattice: 
Mg+^, Nb^'^, and a small number of Pb+"' (with relative 
charges —2, +1, and 0, respectively). The prediction of 
the present model for this case is again just that indicated 
in Fig. 1(c): only a small amount of Pb"*"* (i.e., a small 
value of x) is needed to give rise to the weak 1:1 order 
observed in PMN and PMT. Interestingly, the existence 
of Pb+'* on the B sublattice of ni/2VIi/2 and nii/2Vi/2 
heterovalent alloys, e.g. Pb(Coi/2Wi/2)03 (PCW) and 
PST, will not lead to a change of ordering type. This 
is clearly shown in Figs, l(a-b): whether x is small or 
exactly equal to zero, the 1:1 order is the only order- 
ing predicted by our model for ^(i-x) /2^'^ x'Vl(^i-x) /2 and 
III(i_2:)/2lVxV(i_a.)/2 ternarics. In other words, unlike 
in PMN, the inclusion of a small amount of Pb+^ ions 
in PST may change the strength of the ordering but will 
not change the nature of the order itself. 

While we have chosen to focus above on non-magnetic 
perovskite insulators with B-site substitutional disorder, 
our model is by no means limited to this case. In fact, the 
substitutional "impurity" may even be a vacancy. Denot- 
ing the latter by [ ], our model predicts that the ground 
state of (Cai/2[ ]i/2)Ta03 should adopt the 1:1 rocksalt 
order of the Ca atoms and vacancies on the A sites. This 
is exactly the type of LRO observed experimentally in 
this system [|l9|. The extensions to magnetic perovskite 
insulators or to simultaneous A-site and B-site disorder, 
with possible incorporation of spin contributions or co- 
valent effects into the model, are promising avenues for 
future work. 

In summary, we have demonstrated that a purely ionic 
model is able to reproduce the various types of B-site 
compositional LRO observed in a wide variety of binary 
and ternary perovskite alloys. The remarkable success 
of this approach strongly suggests that the main driv- 
ing mechanism responsible for the LRO occurring in het- 
erovalent perovskite alloys is simply the electrostatic in- 
teraction between the different species. We also raise the 
possibility of a new mechanism in which a small amount 
of Pb"*"* on the B sublattice may be responsible for (or at 
least contribute to) the weak l:l-order X-ray reflections 
observed in some Pb relaxor systems. 
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